Boundary Conditions 


To create a stellar model, one must solve four interrelated differ- 
ential equations, which describe how the radius, pressure, lumi- 
nosity, and temperature change with mass. To do this, we need 
four boundary conditions (to start off each differential). The 
boundary conditions for mass and luminosity are obvious: 


r(m=0)=0 (4.1.1) 


Lim = 0) =0 (4.1.2) 


Unfortunately, we do not have any constraints for the central 
temperature and pressure. For these variables, we must apply 
boundary conditions at the surface of the star. To first order, 
we will define the stellar “surface” as the place where the optical 
depth, 7 = 1. At this location, m = Mr, and we can apply the 
boundary conditions 


P(m = Mr) = 0 (4.1.3) 


T(m= Mr) =0 (4.1.4) 


This approximation isn’t too bad, considering the temperatures 
and pressures in the interior of the star are many order of mag- 
nitude higher. 


In practice, for stars hotter than the Sun, the outer boundary 
condition matters very little, as once inside the star, the tem- 
peratures and pressure quickly converge to similar values. For 
cooler stars, better approximations are sometimes needed, and 
these can make use of relations between luminosity, temperature, 
radius, pressure and opacity. For example, one can use 


Lr = 4r R oT Se (4.1.5) 


(which actually defines effective temperature) and realize that 


the light we see must come from some finite optical depth 
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and note that 


Setting these two equations equal yields 
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P(R) = goT/R = 


where, again, T ~ 1. (Actually, it’s 2/3.) 


(4.1.6) 


Mean Molecular Weight 


The thermodynamic relations between P, p, and T, as well as the 
calculation of stellar opacity requires knowledge of the system’s 
mean molecular weight (defined as the mass per unit mole of 
material, or, alternatively, the mean mass of a particle in Atomic 
Mass Units). Recall that a mole of any substance contains N4 = 
6.02252 x 107° atoms. Thus, the number density of ions is related 
to the mass density, p, by 


n= — = — (5.1.1) 
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where m, is the mass that is equivalent to 1 A.M.U. If the mass 
fraction of species 2 is x;, then its number density is 
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where A; is the atomic weight of the species. The number density 
of all ions in a volume of gas is then 


nt = m= pNad | 
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To compute the contribution of (massless) electrons to the mean 
molecular weight, let Z; be the atomic number of species 7, and f; 
be the species’ ionization fraction, i.e., the fraction of electrons 
of i that are free. The number density of electrons is therefore 


ne = pNay~ (3) Zi (5.1.4) 


or 
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Note that in the case of total ionization (f; = 1), this equation 
simplifies greatly. Since Z;/A; = 1 for hydrogen, and ~ 1/2 for 
everything else, 


u= (x+) =(x+ 95%)" 2 
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(5.1.6) 


From the definitions above, the total number density of particles 
is N 
n = Ng F Ne = A 


where the mean molecular weight is defined as 


p= = + +| (5.1.7) 


The Stellar Opacity 


The mean absorption coefficient, x, is not a constant; it is de- 
pendent on frequency, and is therefore frequently written as «p. 
Inside a star, several different sources of opacity are important, 
and each has its own frequency dependence. To compute the ap- 
propriate “average” opacity, k recall that Fick’s law states that 
the flux at any point is related to the energy density U by 
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and that, when integrated over all energies, 


AacT? dT 
F = ——— — (3.1.4) 
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But how is & related to k? To answer this question, we start by 
explicitly writing in the frequency dependence of U, i.e., 
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and compute the monochromatic flux distribution 


The total flux, integrated over all frequencies is then 
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Comparing (6.1.1) and (3.1.4), we see that 
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50 (6.1.3) 
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The average opacity that you should use is therefore that which 
is harmonically weighted over the temperature derivative of the 
Planck curve. This is called the Rosseland mean opacity. 
Note that the function weights the high frequencies more than 
the low frequencies. To see this, we can take the temperature 
derivative of the Planck function 
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where x = hv/kT, and find its peak by taking the derivative. 
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Mathematically, this works out to x ~ 3.83. In other words, the 
high energy photons count much more than low energy photons. 


Although in real models, the Rosseland mean opacity is com- 
puted numerically, it is instructive to evaluate the form of & 
when ky xv”. From the definition of the Rosseland mean 
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Thus, if the opacity law is kx x v`”, the temperature dependence 
of the Rosseland mean opacity will be k x T7”. This follows 
since the mean opacity must drop as the number of high energy 
photons increases. 


All = 
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ELECTRON SCATTERING 


When free electrons have thermal motions that are much less 
than their rest mass energy, i.e., Mec? >> kT, or T < 5.9 x 
10° K, the cross-section for scattering of photons off electrons is 
frequency independent. The opacity per unit mass of material in 
this case is 
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where ce is the Thomson cross-section of the electron 
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and ne is the density of electrons 


ne = pNay~ (z) fiz (5.1.4) 


If the gas is fully ionized, i.e., if the ionization fractions f; ~ 1 (as 
it will be in stellar interiors), and if we approximate the atomic 
weight (A;) of a metal as twice its atomic number (Z;), then 


ne=PNa et w pNa {T4545} 
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where X, Y, and Z are the mass fractions of hydrogen, helium, 
and metals, respectively. Since X +Y + Z = 1, 


1 1 
Ne = pNasy=+=X (6.1.6) 
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and 


Ke = Naoe(1+ X)/2 = 0.2 (1 + X) cm? g (6.1.7) 


Note again that as long as the ionization is complete (and matter 
is not degenerate), the result is independent of temperature and 
density. 


FREE-FREE ABSORPTION 


The energy emitted per unit volume by a single electron scattered 
by ions of charge Z (as a function of velocity, v) is 
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where e is the charge of the electron, and gg is the gaunt factor. 
To get the emissivity of thermal bremsstrahlung as a function of 
frequency, we integrate this over a Maxwellian distribution, noting 
that to get a photon with frequency v, the incident electron must 
have mev? /2 > hv. Thus, Vmin = (2hv/me)!/? and 
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Next, to go from emission to absorption, we can use the fact that 
in thermodynamic equilibrium, Kirchoff’s law must hold, hence 
the source function is S, = j,/a, = B,(T), where a,, is the 
absorption coefficient per electron (with units of inverse distance). 
Thus 
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or, upon substituting p = Nete/Na and p = nrur/Na 
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The first term in the equation gives the total absorption coefficient; 
the exponential term at the end represents a negative contribution 
from stimulated emission. 


If we combine the constants, letting 


4e® or \ i? 
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then 
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To calculate the Rosseland mean free-free opacity, the inverse of 
this expression, weighted by the derivative of the Planck Function, 
must be integrated over all frequencies. This is tedious, but note 
that if we neglect the effects of stimulated emission, the frequency 
dependence of the opacity is y~°, which, through (6.1.5), implies 
a temperature dependence of T~?. This, plus the T~!/? depen- 
dence which is already in the equation, gives a total temperature 
dependence of T~7/?. Numerically, 
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Opacities of the form k x pT are called Kramers opacities. 


BOUND-FREE ABSORPTION 


For a hydrogenic atom, the cross section for absorption from state 
(n,l) to the continuum is 


SG 2 rimeeZ* gre = 64rtmee!? Z4 grt 
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(6.1.12) 


where the energy of the incident photon must be great enough to 
ionize the atom, t.e., 


From this, it is clear that the opacity law for a single species of 
atom will have a number of absorption “edges” below which the 
absorption cross section will drop quickly. The total bound-free 
opacity, of course, will be the sum over all elements, Z, all ioniza- 
tion states, 7, and and all excitation states, n. In the hydrogenic 
ion approximation, the opacity is 


E Solan, l) (Fe). (saz). = (6.1.13) 
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This summation, along with the process of taking the Rosseland 
mean, will act to smear out many of the edges. 


The bound-free opacity must be computed numerically, but we can 
get a rough idea of how bound-free opacity behaves by making the 
assumption that in a star’s interior, all ions will be very highly 
ionized and thus will have at most one electron. The Saha equation 
states that the number of atoms in ionization state į is related to 
the number in ionization state t + 1 by 


Ni _ Ne ui(T) h? K eXi/kT 
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where u(T) are the partition functions and x; is the ionization 
energy. If the gas is highly ionized, Nj4; ~ 1, ui41 = 1, and 


u; = 2, and 
h2 3/2 JkT 
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(7) 


We can then estimate the number of these atoms in excitation 
state n from the Boltzmann equation. For hydrogenic ions, the 
statistical weight of a level is n°, so if we let Xin = Xi + Xn, then 
the total number of absorbing ions is 
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The total absorption coefficient is therefore 
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and the opacity is 
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where we have dumped the many constants into the variable C 
(and have ignored the stimulated emission term). 


Now examine the factor x;,,/kT. If for a particular n and Z, 
Xin < kT, then that term will contribute little to the opacity, 
since the species will have too few bound electrons. Similarly, if 
Xin > kT, there will not be many photons of the requisite energy 
to be absorbed, and the contribution to the opacity will again be 
small. Thus, only those terms with Xin ~ kT will contribute. 
This leaves kx x pT~'/2v—%, which is the same as free-free opac- 
ity. Bound-free absorption will therefore also follow the Kramers 
opacity law. 


A rough numerical estimate for the Rosseland mean for bound-free 
opacity is 


pew 10? 74x) pr Vom g" (6.1.15) 


Note that this is a factor of ~ 100 greater than the Rosseland 
mean for free-free opacity. 


BOUND-BOUND ABSORPTION 


At extremely high temperatures, T > 10° K, the radiation field will 
consist mostly of ionizing photons, hence bound-bound absorption 
will be minimal (< 10%). However, at cooler temperatures the 
onset of bound-bound transitions in the UV and far UV can double 
the stellar opacity. 


Obviously, the computation of bound-bound opacity must be done 
numerically, since it involves the summation of millions of individ- 
ual absorption lines. It should be noted, however, that in the 
interior of a star, the pressure (and thermal) broadening terms 
will be substantial, so the net result will be continuum absorp- 
tion, rather than a series of absorption lines. Also, the expression 
for the absorption coefficient of hydrogenic bound-free transitions 
(6.1.12) is the same as that for bound-bound transitions. Thus, we 
can expect the same ~ v~° dependence and therefore a Kramer- 
like opacity law. The strength of the total absorption, however, 
will be critically dependent on temperature. 


H- OPACITY 


A free electron passing by a neutral hydrogen atom can produce 
a dipole moment in the atom. At that time, an extra electron can 
attach itself to the system; the result is an H~ ion. This ion is 
fragile — it has an ionization energy of only 0.754 eV, and a large 
cross section for absorption, either through bound-free absorption 


H` +hv © H +e (v) 
or free-free absorption 


H` +e (v)+hv =H +e (v) 


In the envelopes of cool stars, H~ can be the dominant source of 


opacity. 


The contribution of H~ to the star’s opacity will depend on the 
number of H~ ions, which can be calculated from the Saha equa- 


tion. ay 
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where upg- = 1, up, = 2, and y = 0.754 eV. If x is the fraction 
of ionized hydrogen (x ~ 0 for cool stars), and a, is the frequency 
dependency of H~ absorption, then the opacity due to H7 is 
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Note that the total opacity is directly proportional to the electron 
density. In extremely metal poor stars, the free electrons must 
come from ionized hydrogen, hence the opacity will rise with tem- 
perature, until the number of neutrals starts to decline. However, 
for normal metal rich stars, common, low ionization metals (such 
as Ca, Na, K, and Al) contribute to the supply of free electrons, 
thereby lessening this dependence. Consequently, H~ opacity is 
much stronger in metal rich stars than in metal poor stars, and is 
stronger in dwarfs than in giants (due to the pressure term). 


The determination of the frequency dependent cross-section of the 
H- ion involves the quantum mechanical computation of overlap- 
ping wave functions. It is not a power law: it increases linearly at 
short wavelengths until ~ 8500 A, then declines until the ioniza- 
tion threshold at ~ 16,500 A. After that, that opacity increases 
again due to free-free opacity. As a result, the Rosseland mean of 
H~ opacity is not a Kramers law opacity. Very roughly 


Rg- ~ 2.5 x 10731 (4/0102) p1? T? cm? g7! (6.1.17) 
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Bound-free and free-free opacity from H7 at T = 6300°K. Ordinate: 
cross-section { x 104°] per neutral H atom and per unit electron 
pressure p, = AAT; abscissa: 4/1000 where 2 is in A, 


The combined bound-free and free-free opacity from the H~ ion 
at T = 6300 K. The y-axis plots the cross-section (x107~°° per 
neutral H atom per unit electron pressure, Pe = nekT). From 
Mihalas, Stellar Atmospheres (2nd edition), pg. 103. 
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This figure illustrates the importance of the various opacities, as 
a function of temperature and density, for a normal (Population 
I) mix of abundances. The lines are labeled in units of electron 
scattering opacity, Ko = Ke = 0.2(1+ X). From Hayashi, Hoshi, 
& Sugimoto 1962, Suppl. Prog. Theo. Phys. Japan, Vol. 22. 
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The Los Alamos National Labs radiative opacities for X = 0.7, 
Y = 0.28, and Z = 0.02. The dashed line shows the half-ionization 


curve for pure hydrogen. 
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The Rosseland mean opacity (in cm? g~') as a function of density 
(in g cm~?) and temperature for X = 0.739, Y = 0.24, and Z = 
0.24. The range of values plotted is applicable for the outer regions 
of stars; the dots represent values in the solar model. 


